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^ • Abstract 

We consider the two-dimensional quantum field theory of a scalar field self-interacting 
via two periodic terms of frequencies a and (3. Looking at the theory as a perturbed 
Sine-Gordon model, we use Form Factor Perturbation Theory to analyse the evolu- 
tion of the spectrum of particle excitations. We show how, within this formalism, 
the non-locality of the perturbation with respect to the solitons is responsible for 
their confinement in the perturbed theory. The effects of the frequency ratio a/ j3 
being a rational or irrational number and the occurrence of massless flows from 
the gaussian to the Ising fixed point are also discussed. A generalisation of the 
Ashkin- Teller model and the massive Schwinger model are presented as examples 
of application of the formalism. 
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1 Introduction 



For its attractive features, both at classical and quantum level, the (1 + 1) -dimensional 
Sine-Gordon (SG) model has always played a prominent role in the investigation of the- 
oretical non-perturbative aspects of quantum field theory (QFT) [1 — 7] as well as in 
the concrete discussion of innumerable physical systems: equations of Sine-Gordon type 
are in fact relevant in quasi one-dimensional charge density waves of strongly correlated 
electrons and in statistical mechanics aspects of one-dimensional quantum chains 0, in 
self-induced trasparency effects of non-linear optics [10, 11 1 or spin wave propagation in 
the quantum liquid phase of 3 He [12 , just to mention a few. The lagrangian of the model 
may be written as 

C SG = \{d^f + n cos(/^) . (1.1) 

In many of the physical applications, the mapping onto the pure SG model is, however, 
only an approximation and therefore is expected to hold only in some special regions of 
the physical parameter space. A more complete and refined description of several physi- 
cal effects often requires that a lagrangian with two or more periodic interaction terms is 
considered. A simple example is provided by a system of two Ising models coupled both 
by thermal and magnetic operators, system which may be considered as a generalization 
of the Ashkin- Teller modelQ. Other examples, arising from different physical contexts, 
are given by the aforementioned cases of spin waves propagating in anisotropic magnetic 
liquids or ultra-short optical pulses propagating in resonant degenerate medium |T3) . The 
purpose of this paper is to analyse the quantum field theory aspects of this kind of gen- 
eralisation of the SG model and in particular the case in which two periodic interactions 
are present, i.e. to study the QFT defined by the lagrangian^] 

C = -{d^(p) 2 + fi cos j3ip + A cos(a<£> + 5) . (1.2) 

Apart from the practical motivations we have just mentioned, there are also some the- 
oretical issues which make the analysis of (|1.2j ) an interesting problem. Let us try to 
describe them. 

The addition of a new interaction term to the lagrangian fll.l] ) spoils one of the crucial 
properties of the SG model, i.e. its integrability. As a consequence, the resulting two- 
frequency model ( |1.2j ) is not exactly solvable and perturbation theory is basically the 
only available tool for a quantitative study. Of all the different perturbative approaches 
that could be applied to the problem at hand, the one where the new periodic term is 
viewed as perturbing a pure SG model is probably the least conventional but the most 



1 This model will be discussed in some details in Section 5 of this paper. 

2 Meaning and restrictions on the parameters entering eq. (1.2) will be the subject of Section 2. 
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convenient. Perturbation of integrable models was discussed in Ref. [15]: its theory relies 
on the fact that the matrix elements of local operators between asymptotic states of the 
unperturbed (integrable) theory are exactly computable. The corrections to the masses of 
the particles of the unperturbed theory and to the scattering amplitudes can be computed 
in terms of the matrix elements of the perturbing operator in much the same way as in 
quantum mechanics. Since the matrix elements of the local operators are usually known 
as Form Factors, we refer to this perturbative approach as Form Factor Perturbation 
Theory (FFPT). The Form Factors of the SG model were computed in Refs. |5], ^ 
and consequently the technical information required for these perturbative calculations 
is available. However, our interest is not to perform specific perturbative computations, 
but rather to illustrate some general features of the perturbed model and the use of the 
FFPT to describe it. 

We are particularly interested in the spectrum of the physical excitations. The invari- 
ance of the Lagrangian ( |1 . 1| ) under the discrete transformation (p — > (p + 27rg//3, q G Z 
implies that the elementary excitations of the pure SG model are solitons interpolationg 
between adjacent degenerate minima of the potential. Extra interactions will generally 
remove the degeneracy among the different classical vacua and therefore these topologi- 
cal excitations are expected to disappear from the spectrum of asymptotic states of the 
perturbed theory. Therefore it is essential to ascertain whether the FFPT accounts for 
such a dramatic evolution in the particle spectrum of the theory. As we will show, the 
answer is indeed affirmative and the phenomenon of confinement of the original topologi- 
cal excitations emerges in the FFPT as a consequence of the mutual non-locality between 
the perturbing operator and the field which interpolates the confined particle. Notably, 
the non-local mechanism responsible for the confinement of the topological excitations 
proves to be a general feature of many two-dimensional theories: it causes, for instance, 
the disappearence of the kinks from the spectrum of the low-temperature phase of the 
Ising model as soon as an external magnetic field is turned on [15], |I6 |. 

The role of the quantum vacuum of the theory is particularly relevant for the analysis 
below. In the pure SG model, the quantum vacuum is selected from the infinite number 
of classical equivalent vacua through the spontaneous symmetry breaking mechanism. At 
the operational level, this involves selecting arbitrarily one of the classical minima as the 
quantum vacuum and constructing the space of asymptotic states as excitations above it. 
Physical quantities such as scattering amplitudes, Form Factors, etc. are then computed 
as matrix elements on this space of states. In the perturbed theory, the degeneracy of 
the classical vacua is explicitely broken and the quantum vacuum is uniquely defined. A 
fundamental requirement for the application of the FFPT is to ensure that, when the 
perturbation is adiabatically switched off, the vacuum matches the one previously chosen 
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for the computations in the unperturbed theory. This "adiabatic condition" corresponds 
to the general requirement that the vacuum of a quantum field theory in the presence 
of a spontaneously broken symmetry must be identified by taking the limit in which a 
symmetry breaking perturbation is made to vanish. 

The vacuum of the perturbed theory is however unique and the above considerations 
apply, only if the periods of the two interaction terms are incommensurate. In the case 
of rational frequency ratio, the survival of a subset of degenerate classical vacua in the 
perturbed potential must be given appropriate consideration. It is already clear from this 
observation that the theory will exhibit an extremely irregular behaviour as a function 
of the frequency ratio. 

In the non-integrable model ( |1.2| ), the two periodic interactions clearly play a sym- 
metric role and each of the cosine interactions can be viewed as a deformation of the 
integrable theory defined by the other: as discussed in more detail in Section 3, there is 
a dimensionless variable 



which characterises the two perturbative regimes of the model ( |1.2|) . The first regime is 
obtained in the limit 77 — > 0, while the other is reached for rj — > 00, the latter by just 
swapping the role played by the two operators. The possibility to control the changes 
of the spectrum in both perturbative limits allows one to deduce interesting information 
about its evolution in the intermediate, non-perturbative region. In certain cases, there 
may be a topological excitation in one limit which is no longer present in the other. When 
this is the case, the very nature of topological excitations requires that a change in the 
vacuum structure of the theory takes place somewhere in the non-perturbative region, 
namely that a phase transition occurs^. Lines of phase transition are then expected to 
appear in the multi-frequency SG model for particular values of the parameters: they 
generally correspond to Renormalisation Group (RG) trajectories along which the system 
flows from the gaussian fixed point with central charge C = 1 to the Ising fixed point 
with central charge C = 1/2. 

The paper is organised as follows. The next section reviews several results of the pure 
SG model which will be used later. The two-frequency model is introduced in Section 3: 
some restrictions on the parameters of the theory dictated by conformal perturbation 
theory and the aforementioned adiabatic condition will be discussed. The Form Factor 
Perturbation Theory, the relation between non- locality and soliton confinement, the evo- 
lution of the particle spectrum and the occurrence of phase transition are analysed in 

3 This phenomenon was discussed in these terms long ago by S. Coleman in his classic study of the 
massive Schwinger model |2q| . 



= ^ -(1-A a )/(1-A p ) = A ^-(8^ Q 2 )/(8^-/? 2 ) 
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Section 4. The system of two-layer Ising model, alias the generalized Ashkin- Teller model, 
and the massive Schwinger model are presented as examples of specific applications of the 
formalism in Section 5. Summary of the results and conclusions are in Section 6. The pa- 
per has also two appendices: the first contains a derivation of the soliton-antisoliton form 
factors of exponential operators at the reflectionless points of SG model while the second 
presents an analysis of soliton confinement within the semiclassical approximation. 



2 Few Notions of the Sine— Gordon Model 

In this section we briefly review some basic features of the pure SG model and establish 
several notations which will be used in the following. The SG model is a relativistic 
theory in (1 + 1) dimensions for a scalar bosonic field <£>(x) with the lagrangian given in 
eq. fll.lj) . At the classical level, the model possesses an infinite degeneracy of the vacuum 
states which can be labelled by a relative integer q: finite energy configurations must 
necessarely approach one of the vacuum values ^S-q- for x — > — oo and another one 2 ^-q+ 
for x — > +oo. The difference Q = q + — q_ plays a role of topological charge which 
characterises inequivalent sectors of the theory. As well known, the model is completly 
solvable both at the classical and the quantum level. At the classical level, all solutions 
of the equation of motion can be computed by means of the inverse scattering method, as 



reviewed for instance in [14] . Building blocks for expressing all the classical solutions are 
the static soliton and antisoliton configurations with Q — ±1: their explicit expressions 
are given by 

4 

ip s ,s(x) = -arctan[exp(±/3 v //Ia;)] , (2.1) 

where + refers to the soliton and — to the antisoliton. Multi-solitons are then obtained 
by the non-linear superposition principle provided by the Backlund transformations. For 
the purposes of this paper, it may be convenient to think of soliton configurations as step 
functions with jumps ±27r//3 (Figure 1). 

At the quantum level, a normal ordering of the cosine interaction term in the la- 
grangian ( |1 . 1| ) is sufficient to get rid of the only diverging diagrams (tadpoles) in the 
theory. From the conformal perturbation theory viewpoint, the SG model can be consid- 
ered as a deformation of the gaussian fixed point action 

Acaussian = J d 2 X^{8 v (ff , (2.2) 

(with central charge C — 1) by the operator cos (3(p with conformal dimension = 
j3 2 /8n. The latter is relevant for f3 2 < 8ir, and in this range the theory presents a massive 
phaseQ The exact particle spectrum of the theory |fj, || consists of a pair of soliton and 
4 The limitation (3 2 < 8tt will be always assumed to hold in the following. 
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antisoliton excitations with equal mass M and a number N 
bound states B n (breathers) with masses 



of soliton-antisoliton 



m k = 2Msm[^Pj , k = 1, 2, . . . , < - , (2.3) 



where 



e = ^ • (2.4) 



7l(3 2 

The elementary field tp entering the lagrangian ( |1 . 1| ) may be regarded as associated to 
the lightest breather B\. Note, however, that the fc-th bound state crosses the soliton- 
antisoliton threshold when £ = n/k, and for £ > n/k it leaves the physical spectrum. 
Hence, for £ > tt, the spectrum consists of soliton and antisoliton only and does not 
contain any particle associated to the field (p. 

An additional fundamental result about SG model is its exact equivalence with the 
massive Thirring model [|l|], namely with the theory of a self- interacting Dirac fermion 
defined by the Lagrangian 

Cmtm = ^i-fdp - m)^ - |(^) 2 • (2.5) 

The correspondence between the two models is specified by the bosonisation rules and 
by the following relation among the couplings of the two models 

T 2 

The exact S"-matrix relative to factorizable and elastic scattering processes of the SG 
model has been determined in 0. We refer the reader to this paper for all details about 
this subject. The matrix elements of local operators between the asymptotic states of 
the SG model have also been computed by using the so called Form Factor bootstrap 
approach |5], |6| . Recently, Lukyanov used a technique based on free field representation to 
compute the Form Factors of the exponential operators e %aip [^]. Since these Form Factors 
are of particular interest in the present paper, we show in Appendix A how they can be 
easily computed in the ordinary bootstrap framework, at least for the special values of 
the coupling £ = n/k where the S'-matrix is reflectionless and the technical problem gets 
considerably simplified. 



9 = • (2.6) 



3 The Two-frequency Sine-Gordon Model 



In this section we will derive and illustrate some limitations on the real parameters a, (3 
and S entering the quantum action 

"1 



A 



d x 



^(d^ip) + cos f3ip + \cos(aip + 5) 



(3.1) 
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We can always assume that a, /3, \i and A are positive, since this situation can be arranged 
by using appropriate redefinitions of ip and 5. 

There are at least two ways of viewing the above theory, both of which provide some 
useful insight to its physical content. The first regards (|3.1|) as a perturbation of the gaus- 
sian action ( |2.2| ) by means of the two scaling operators 0\ = cos /3(p and O2 = cos (atp +5), 
with conformal dimensions = f3 2 /8n and A a = a 2 /8n, respectively. The coupling 
constants /x and A have scaling dimensions [mass) 2<yl ~ Afi ^ and [mass) 2<yl ~ Aa \ respectively 
and therefore the parameter 77 which was defined in ( |1.3| ) forms a dimensionless combina- 



tion which labels the different RG flows originating from the gaussian fixed point. The 
relevant nature of both operators then implies 

^ < 5 (3-2) 
(3 2 < 8tt . V ; 

An additional constraint on the two parameters a and (3 derives from the renormalizability 
of the action (|3.1|) , an aspect which can be efficiently studied by conformal perturbation 
theory [17]. To illustrate this point, let us briefly discuss the pattern of renormalisation 



of a scaling operator in a theory obtained by perturbing a conformal actionf] by a 
relevant operator <fi of conformal dimension A 

A g = AcFT + g J d 2 x<P(x) . (3.3) 

Let X denote a generic product of operators and consider the usual perturbative expan- 
sion of the correlator 

(X$(0)> = (X$(0)) C FT + g [ d 2 x(X<b(Q)<P(x)) C FT + 0{g 2 ) . (3.4) 

e and R are ultraviolet and infrared cutoffs, respectively, and the correlators on the right 
hand side of the above equation are computed at the conformal point. The integral in 
Eq. (|3.4|) is ultraviolet divergent only if the conformal OPE 

$ (°) = E^i x i 2{Afc_A ^ A) ^(°) ( 3 - 5 ) 



[X 



contains operators with conformal dimension A*, such that 

7/c = A fc - A$ - A + 1 < . (3.6) 

If this is the case, a first order UV finite correlator (X $'(0)) is obtained by defining the 
renormalised operator 

$' = <5> + gY,he 2 ^A k + 0(g 2 ) , (3.7) 



5 We consider for simplicity the case of perturbation by a single scaling operator, but the following 
considerations are easily extended to more general situations. 
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where bk = —irCi^/jk an d the sum runs only over the operators A k satisfying the condi- 
tion ( |3.6| ). The last equation shows that the renormalisation of UV divergences induces a 
mixing of the operator with a finite number of operators having smaller conformal 
dimensions. 

The above general considerations are easily applicable to perturbations of the gaussian 
action, with the spectrum of primary operators spanned by the exponentials e iCttp ( x ' of 
conformal dimension A a = a 2 /87r and their OPE given by 

e iai<p(x) e iot2<p(0) _ | a ,jaia2/27r e »(ai +03)^(0) (3 8) 



Using eq. (|3.6| ), it is easy to see that the two perturbing terms in ( [3.1|) do not mix with 
other operators (so that no extra terms are generated in the action fl3.1| )) if the parameters 
a and (5 satisfy the inequality 

ap<4n . (3.9) 



If the above condition is not fulfilled, the action (|3.1|) results unstable under renormaliza- 
tion, i.e. additional counterterms e ±l< - Q_ ^ should already be included at the first order 
in g so as to have a consistent theory^. 

Let us now consider the second approach to the QFT ( |3.1| ) - the one we will mostly 
adopt in this paper - which consists of viewing it as a deformation of a pure Sine- 
Gordon action. Although the two interaction terms play a completely symmetric role, in 
the following we will assume A cos(a(p + 5) to be the perturbing term and derive some 
restrictions on the parameter 5. 

When A = 0, the pure SG potential has a periodicity 27r / /3 and therefore exhibits an 
infinite number of classical degenerate vacua placed at (p — 2nq//3. At the quantum level, 
however, spontaneous symmetry breaking selects a unique quantum vacuum, which by 
convention we choose to be at the origin, (p = 0. Let us now switch the perturbation on. 
Our analysis begins by considering the case when the frequency ratio 

w = | (3-10) 

is a rational number, m/n, with m and n coprime positive integers. In this case, the new 
potential acquires a periodicity 2im/ (3. For A very small, the perturbed potential presents, 
within the first period, n relative minima at (fk = j3+0(X), k — 0,1, ... , n—1 (Figure 
2). Let the absolute minimum of the perturbed potential be located nearby <p~ k = 2nk/{3. 
According to the adiabatic condition, we have then to shift 

tp{x) — >ip-2nk/(3, (3.11) 



6 In the particular case a = f3 this simply leads to the introduction of an inessential additive constant 
in the action (ti.lf). 
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in order to move the vacuum back near to the origin: the action (|3.1|) is intended to 
precisely address to this situation. Hence, the only values of 5 which lead to inequivalent 
quantum field theories - those which do not differ by a simple shift in the location of the 
vacuum - are those given by 

\S\ < I . (3.12) 

n 

These are in fact the values for which the absolute minimum of the perturbed potential 
remains located close to tp = 0. At the perturbative level in A, eq. ( |3.12| ) is simply derived 



from the requirement that the variation of the potential at the minimum close to tp — 
is larger (in absolute value) than at the other (n — 1) minima. When the bound ( |3.12| ) is 
saturated, one of the (n — 1) minima and the one close to the origin degenerate (Figure 
3). 

Let us now consider the case in which the ratio u of the two frequencies is an irrational 
number. Since any irrational number can be approximated with arbitrary precision by 
the sequence of rational approximants provided by its continued fraction expansion, this 
case may be viewed as a particular limit of the previous situation. Let 

m,- 1 . 

(3-13) 



Ui ri + 



<■_>- 



be the z-th approximation of the irrational number uj. Both rrii and n,i become larger 
by increasing the index i of the sequence^. Since the restriction (|3.12|) only involves 
the denominator n« of the rational number, the range of values of 5 which give rise to 
inequivalent theories shrinks consequently to zero in the irrational case. 

In conclusion, the family of RG trajectories in the A/x-plane which are labelled by 
the dimensionless parameter rj acquires, somehow, a "fractal" nature as a function of the 
frequency ratio u: while no further parametric dependence is present for irrational values 
of u, on the contrary a multiplicity of trajectories - labelled by the inequivalent values 
of 5 ( p.!2|) - results for each rational value of the frequency ratio. 



The varied behaviour of the theory ( |3.1| ) - depending on whether w is a rational or 
an irrational number - can be understood in a number of ways. In the presence of an 
irrational frequency ratio, the perturbed potential has a very irregular shape: it is no 
longer periodic and possesses an infinite number of inequivalent relative minima. The 
values of the potential at these minima form a continuous spectrum above the absolute 
minimum value ^nin = ~~ A 4 — ^- Suppose that the absolute minimum lies at tp — tp' for 
a given value of 5 (e.g. tp' = for 5 = 0). An infinitesimal variation of the parameter 5 
induces a corresponding variation of the potential at tp' . Since the values of the potential 



7 A well known example of this statement is provided by the sequence Fi/Fi+i of the Fibonacci 
numbers which converges to the golden ratio (vo — l)/2. 
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at the minima form a continuous spectrum, the absolute minimum will no longer be 
located near ip' but will jump to another point ip" which is generally far from <p'. Hence, 
a variation of 8 simply amounts to a shift of the vacuum from ip' to tp". This in turn 
means that 8 is an irrelevant parameter in the case of irrational frequency ratio and can 
be safely set to zero. 

The above considerations can also be rephrased in more geometrical terms, as follows. 
Let us introduce an additional field and write the potential of the lagrangian ( |1.2[) 

as 

V((p,Q = -fJ- cos (3<p- A cos(«C + S) . (3.14) 
The original theory is of course recovered by the identification 

( = <p . (3.15) 

Independently from the nature of the ratio u, the potential V(ip, () presents a double 
periodicity (^-, as a function of the two variables p and £. Hence, this function 

is naturally defined on a compact torus T — C) : ~ | < V 9 < f > — ^ < C < f }> an d 
in this domain has an absolute minimum located at P* = (0, — -). To this point, the 
rational or irrational nature of u has played no role. Only by enforcing the identification 
( |3.15| ) does the difference between the two cases emerge. Note in fact that on the torus T 
eq. (|3.15| ) defines a family of straight lines which all have slope 1 (Figure 4). These straight 
lines may be viewed as the motion of a geometrical point bouncing on the boundaries of 
the torus and reappearing on the opposite sides via its periodicity. When uj is a rational 
number m/n, the motion of this point repeats over along the finite number n + m — 1 
of these diagonal lines (Figure 4. a). The motion may pass through the minimum P* or 
remain at a finite distance from it, which depends on the value of 8. It is however clear 
that in the rational case all values \8\ < - are inequivalent and that when 8 exceeds this 
bound we can go back to the case already considered by simply redefining the funda- 
mental domain T of the torus. In the irrational case, on the contrary, the motion of the 
geometrical point is ergodic and therefore spans the entire torus T (Figure 4.b). Hence, 
it passes arbitrarly close to the minimum P*, regardless on the value of 8. Consequently, 
the theory with 8 ^ is physically indistinguishable from the one with 8 = 0, the latter 
value only reflects the choice of cp = as the vacuum state. 



4 Non-locality, Soliton Confinement and Phase Tran- 
sitions 

In this section we will consider the quantum theory (|3.1| ) as a perturbed Sine-Gordon 
model with the purpose of investigating some features of the evolution of the particle 
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spectrum of the theory as a function of the parameter 7] defined in ( |1.3|) . 



4.1 Form Factor Perturbation Theory and Non— locality 



Perturbation theory around integrable models was discussed in Ref. [TJ|, where it was 
shown how the corrections to the mass spectrum and scattering amplitudes can be ex- 
pressed in terms of the exactly computable Form Factors of the perturbing operator on 
the asymptotic states of the unperturbed theory. For the purposes of this paper we only 



need to point out some of their properties and few basic results of Ref. [[HJ . 
Let be the scalar operator which perturbs an integrable action, 

A = Ant + A J d 2 x^{x) . (4.1) 

One of the first effects of moving away from integrability is a change in the spectrum of 
the theory: the first order correction to the mass of a particle a belonging to the spectrum 
of the unperturbed theory is given by 

5ml ~ 2AF*(*7r) , (4.2) 

where the particle-antiparticle Form Factor of the operator defined as the matrix 

element]^ 

Fl{6) = (OMO)\<0iH0 2 )) , (4-3) 

is introduced. 

In integrable theories Form Factors of a generic scalar operator 0{x) satisfy man- 
ageable functional equations in virtue of the simple form assumed by the unitarity and 
crossing symmetry equations || |6|. For the two-particle case, we have 

Fg(e) = S*(0)Fg(-0) , (4.4) 

Fg(9 + 2m) = e- 2 ^Fg{-6) . (4.5) 

The first of these equations expresses the fact that in an integrable theory the two- 
particle threshold is the only unitarity branch point in the plane of the Mandelstam 
variable s = (p a + Pa) 2 = 4m^ cosh 2 (#/2), the discontinuity across the cut being deter- 
mined by the two-body scattering amplitude S^ b a . 

In the second equation the explicit phase factor e _2l7r7o > a is inserted to take into 
account a possible semi-locality of the operator which interpolates the particle a (i.e. 

8 The standard parameterisation of two-dimensional on mass-shell momenta in terms of rapidities is 
adopted: = m a cosh0i, p\ = m a sinh#i, 6 = Q\ — 62- 
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any operator <p a such that (0|<^ a |a) 7^ 0) with respect to the operator 0(x) f\. We must 
bear in mind that two operators 0\ and O2 are said to be mutually non-local if the 
euclidean correlator (. . . Oi(xi)02(x2) ■ ■ .) is not a single valued function of X\ and x<i- 
In the particular case in which the correlator simply acquires a phase factor under the 
analytic continuation bringing x<i to its original position along a path encircling xi, then 
the operators 0\ and 2 are said to be semi-local with respect to each other. Equation 
( j4.5|) can be regarded as a momentum space version of this statement for the operators 
0(x) and ^„(x). When 7^ = 0, there is no crossing symmetric counterpart to the 
unitarity cut but when jo,a 7^ 0, there is a non-locality discontinuity in the s-plane with 
s = as branch point. In the rapidity parameterisation there is however no cut because 
the different Riemann sheets of the s-plane are mapped onto different sections of the 
0-plane; the branch point s = is mapped onto the points 9 = ±m which become the 
locations of simple annihilation poles, whose residues are given by [Q (see also |18| ) 

- iRes e=±i «Fg{e) = (1 - e* 2i ^)(0\O\0) . (4.6) 

With the above information, let us now proceed to the perturbative analysis of the 
action (|3.1|) . As before, in the following we will refer to the case 77 <C 1, so that the role 
of the perturbing operator \l/ in ( |4.2| ) is played by cos(ap + 5). It is understood that a 
similar analysis can be repeated in the opposite limit rj — > 00, with /i cos flip assumed 
as the perturbation. 



4.2 Fate of the solitons 

The particle spectrum of the unperturbed theory at rj = consists of a soliton-antisoliton 
pair and - for (3 2 < Att - a certain number of breathers. 

Breathers are mutually local with respect to the elementary Sine-Gordon field (p, and 
consequently to the perturbing operator cos(a<^ + 5). Thus, formula ( f4.2p can be safely 
used to explicitely compute the first order correction to the unperturbed masses ( p.^ )! 10 !- 

The situation is quite different for the solitonic sector of the Hilbert space. The 
soliton is an elementary excitation of the unperturbed theory which interpolates between 
two constant configurations of the field tp differing by 2tt/(3. This means that ip — > 
p + 2n/(5 across a soliton configuration. As a consequence, the exponential operator e iaip 
is semi- local with respect to the soliton with semi-locality index 7^ )S = a/ (3. According 
to equation (^4.6| ), the soliton-antisoliton Form Factor of the perturbing operator \& = 
cos{ap + 5) contains annihilation poles in 9 = ±iir whose residues are immediately 
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Consistency of eq. (4^) requires 70, a = — lo,a- 

10r 



3 The Form Factors which are actually needed can be found in (?], |). 
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computed as^ 

- iRes e =±inF? s (6) = [cos 5 - cos(<5 q= 2ira//3)] (0\e iav \0) . (4.7) 

The presence of these poles has a drastic consequence on the spectrum of the solitonic 
sector of the perturbed theory. Indeed, for generic values of a and 5, the non-zero residue 
Q4.71) implies an infinite correction to the soliton and antisoliton masses. This divergence 
may be viewed as the technical manifestation of the fact that soliton and antisoliton 
no longer survive as asymptotic particles of the perturbed theory. Of course, this result 
could be anticipated by observing that for generic values of a, the lagrangian (|1.2[) loses its 
original 27r//3-periodicity, so that solitons and antisolitons become unstable excitations. 
As discussed in more detail in Appendix B, the lifting of the degeneracy of the original 
minima induces a linear attractive potential between the soliton-antisoliton pairs which 
causes their collapse into a string of bound states. 

Let us consider more closely the situation in which uo = m/n, with m and n coprime 
integers, so that the lagrangian dl.2|) retains a 27m//3-periodicity. As just argued, as 
soon as the perturbation A cos(a<£> + 5) is switched on, the solitons of the original theory 
become confined into states with zero topological charge^ but in this case "packets" 
formed by n of the original solitons (or antisolitons) survive as stable excitations. These 
are nothing else but the topological excitations which interpolate between the degenerate 
minima of the perturbed potential^ (see Figure l.c and 2.b). Clearly, a complete space 
of topological excitations for the perturbed theory can be constructed in terms of the 
"n-soliton" and "n-antisoliton" states. 

The above conclusion holds for generic values of 5. However, for specific values of 
this parameter, the single-soliton states of the unperturbed theory may remain stable 
and unconfined even in the presence of the perturbation. In the FFPT framework, this 
corresponds to the case in which the correction to the soliton mass is finite, namely when 
the residue on the annihilation pole ( [4.7D vanishes. Taking into account the condition 
( |3.12| ) satisfied by 5 in the rational case, the residue ( f4.7| ) vanishes when 

\5\ = it/n (4.8) 

provided m and n satisfy the condition 

\kn — m\ — 1 for some k e Z . (4-9) 

11 This result refers to the choice (0 1 1 0) = for the unperturbed theory. In this case the unperturbed 
theory is invariant under the reflection ip — > — <p and the equality (0|e IQ¥, |0) = (0|e -za¥> |0) follows. 
12 This is a statement which holds orn^l and generic values of S. 

13 Such packets, considered as a whole, are local with respect to the exponential operator e la(p since 
their semi-locality phase is e ± 247rnQ // 3 = 
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These are precisely the values of lu and 5 for which the minimum at p = of the 
unperturbed potential becomes degenerate with one of its adjacent minima when the 
perturbation is switched on (Figure 3.b), and therefore the single-soliton excitations 
interpolating between the two minima remain stable. Note, however that it is generally 
impossible to simultaneously cancel both the annihilation poles at 9 = m and 9 = —m\ 
it can be argued that this corresponds to the fact that only one of the adjacent minima 
(the one located at <p — 2ir/j3 + 0{rf) or the one located at <p = —2n/(3 + 0{rj)) can 
become degenerate with the absolute minimun around p = (order rj) . Suppose that the 
degeneracy is realised between the central minimum and its right neighbour. The presence 
of the degeneracy implies that, for the values of the parameters we are considering, the 
n-soliton discussed above actually breaks into the sequence of two asymptotically stable 
excitations: a kink K\ of mass M interpolating between p = and p = 2tc//3, and a kink 
K n -i of mass (n — 1)M interpolating between (p = 2tt//3 and p = 2im/ (3 (here and below 
0(r]) corrections to the masses and positions of the minima are understood). Obviously, 
this pattern is repeated with periodicity 27in/(3. 

More generally, it is not difficult to see that under the action of the perturbation, for 
6 = ETi/n (e = ±1), and for any pair m,n (m ^ 1), the absolute minimum close to p = 
remains degenerate with the minimun located near to tp — 2Tcj//3, with j given by 

j = e{l,2,...,n-l} . 4.10 

m 

There is always a unique integer k for which the above equation is satisfied. In this case, 
the n-soliton breaks into a kink Kj with mass jM and a kink K n _j with mass (n — j)M. 

The denomination of kinks we gave to these excitations requires some clarification. 
Generally, kinks are elementary excitations which interpolate between different vacua 
of a theory with spontaneously broken symmetry. When the structure of the vacua 
is sufficiently non-trivial, suitable ordering prescriptions must be assigned to construct 
acceptable multi-kink configurations. In this circumstance, kinks drastically differ from 
ordinary particles, in terms of which the space of multi-particle states can be constructed 
without any restriction. The situation we are dealing with is not of this kind. In fact, 
due to the simple vacuum structure (alternation of two different kinds of minima), any 
string of Kj and K n -j is an allowed multi-kink configuration. Differently stated, for 
any sequence of Kj and K n -j there exists an unique interpolation pattern among the 
degenerate vacua of the theory^. 

14 Since we are dealing with excitations above the real vacuum of the theory, the interpolation starts 
from ip = 0. Moreover, the sequences KjKj and K n -jK n -j can only correspond to going back and forth 
between a couple of consecutive minima. Then, as long as the alternating minima remain inequivalcnt, 
there is no need to define antikink states. 
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It should be clear that the above considerations regarding the evolution of the particle 
spectrum of the theory ( |3.1| ) have a validity which goes beyond the perturbative frame- 
work in which they have been derived: the stable or instable character of an excitation 
under the action of the perturbation is a qualitative feature which can be extended away 
from the perturbative domain. For example, had we performed our perturbative analysis 
in the limit rj — ► oo rather than in the limit 77 — > 0, we would have established the 
existence of m-solitons rather than n-solitons. These are different names for an unique 
excitation which exists and is stable for any value of 77. Simply, it is intuitively helpful to 
see this object as the bound state of n solitons of the 77 = theory, in one perturbative 
limit, or m solitons of the 77 = 00 theory, in the other. 



4.3 Phase Transition 

Let us apply the above considerations to the case 5 = ±7r/n in which the n-solitons 
"decay" into the sequence of kinks Kj and K n _j. By increasing 77, the masses of these 
particles evolve preserving though their stability, and for 77 very large they can be rein- 
terpreted perturbatively as K\ and K m _i kinks, with I determined by eq. Q4.10D but with 
m and n interchanged. 

An exception to the above smooth evolution pattern arises when m = 1: in this case, 
the n-soliton - which may be seen as "sum" of two kinks K\ and K n _\ - evolves for 
77 — > 00 into a "1-soliton". The latter, up to a mass correction, is simply the soliton 
of the 77 = 00 theory and therefore is completely stable (in this limit the perturbation 
is local with respect to all the excitations of the unperturbed theory for m — 1). Thus, 
a composite topological excitation has surprisingly transmuted into an elementary one 
along the way! 

Such a transmutation necessarily requires the existence of an intermediate critical 
value t] c at which a phase transition takes placep 5 ]. The phase transition presents the 
following pattern. For 77 < r/ c , the n-soliton consists of the sequence of the two kinks K\ 
and K n -\. When 77 — > rj c from below, the intermediate minimum which is in common 
between the two kinks approaches the one initially located at </? = (Figure 5. a). The 
two minima become coincident at 77 = r? c (Figure 5.b) and, at this value, K\ disappears 
as a topological excitation. In particular, its mass (which is proportional to the potential 
barrier between the minima) shrinks from M at rj = to at 77 = r] c . For rj > rj c 
the original composite topological excitation consists of a single massive kink, which 
interpolates between the vacua of the 77 = 00 theory (Figure 5.c). However, the dynamical 
degree of freedom carried by K\ do not disappear across the transition line: it is in fact 
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One could imagine several critical thresholds. The simplest possibility is assumed here. 
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transferred to a particle with zero topological charge which is precisely massless at 77 = 7] c 
and becomes a breather at 77 = 00Q. This tranformation of a topological excitation into a 
non-topological one should not be surprising: we already remarked that, from the point 
of view of the construction of the space of asymptotic quantum states, K\ is actually 
indistinguishable from an ordinary particle. 

The presence of a massless particle in the spectrum of the 77 = r) c theory amounts to 
say that the correlation length is infinite along the RG trajectory labelled by ?7 C . The 
correlation length cannot decrease in the RG flow towards the infrared limit, therefore 
such a trajectory must end into a new fixed point. Since the theory ( [3.1|) is unitary, 
Zamolodchikov's C-theorem Jl9| ensures that such a fixed point is described by a CFT 
with central charge smaller than that of the ultraviolet fixed point, namely C — 1. On 



the other hand, it is known from Ref. |20| that the only unitary CFTs with C < 1 are 
those belonging to the unitary minimal series, with central charges C p — 1 — Q/[p(p+ 1)], 
p — 3,4, . . .. For 7] slightly below i] c , the mass of K\ is much smaller than the mass of 
the other particles of the theory. Thus, the low energy limit is effectively described only 
by the kink K\ interpolating among the vacuum and its adjacent degenerate minimum. 
Stated in a Landau- Ginzburg language, the action ( |3.1| ) can be effectively truncated to a 
spontaneously broken f 4 theory, which is well known to describe the universality class of 
the Ising model. In fact, the one we described above is exactly the Ising phase transition 
with the parameter rj playing the role of the temperature. Then we conclude that the 
RG trajectory labelled by r] c interpolates between the gaussian fixed point with C = 1 
and the Ising fixed point with C = C3 = 1/2 (Figure 6). A model which provides an 
explicit realization of this flow is discussed in the next Section. 



5 Two examples 

This section is intended to illustrate the wide applicability of the ideas discussed in this 
paper through a couple of examples referring to quite different physical situations. The 
first example is taken from classical statistical mechanics and deals with a deformation of 
the Ashkin- Teller model. The second one comes from quantum field theory and concerns 
the massive Schwinger model. Our interest is in exhibiting those features of the two 
models which are relevant for the discussion developed in this paper. The reader is 



referred to the original literature for comprehensive studies on both subjects (see Refs. [^T 
22j for Ashkin- Teller and ^5j for the massive Schwinger model). 



16 Note that the existence of at least one breather in this limit is ensured, for m — 1, by the condition 
a/3 < 4-7T we derived in the previous section. 
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5.1 The Generalized Ashkin- Teller Model 



The Ashkin- Teller model describes two planar Ising models interacting through a local 
four spin interaction. It is defined by the lattice Hamiltonian 

Hat = ]T[J(aX + o\a\) + Ka\a\a\a{\ , (5.1) 

(id) 

where o\ 2 = ±1, the sum is over nearest neighbours and the same coupling J has been 
chosen for the two Ising models (isotropic case). 

A scaling limit description of the model can be formulated in terms of the euclidean 
action of a two-layer Ising system (Figure 7) 

Aat = ^. + ^. +T j A(ei(l) +S2ix))+p J d ^ M) , ( , 2) 

where Jl 1 .^ 9 and £i(x) denote the fixed point action and the energy operator of the z-th 
Ising model, respectively. In this language, the scaling limit of the Ashkin- Teller model 
appears as a CFT with central charge C = 1/2 + 1/2 = 1 perturbed by the two operators 
£ = Ei + e 2 and e = eiE 2 - 

Let us keep at first the two Ising models at their critical temperature (r = 0). Since 
the operator e is marginal (A £ = 1/2 + 1/2), its addition to the action of the two 
critical Ising models does not spoil criticality, but leads to a line of C = 1 fixed points 



parameterised by the values of the coupling p. It is known [21], |22] that, while the spin 
fields <7i and 02 retain along the critical line the conformal dimensions which they have 
at the decoupling point p = (i.e. A CT = 1/16), the dimensionalities of the total energy 
operator £ = E\ + 62 and of the so called "polarisation" operator P = a\(J2 become 
instead p-dependent. Their ratio, however, remains the same as at the decoupling point 

Ap(p) = MO) = 1/16 + 1/16 = 1 

A £ (p) Ae(0) 1/2 4 • 1 ' ' 

Thus, the action fl5.2|) with r / can be seen as a C = 1 CFT perturbed by the operator 
£ with p-dependent conformal dimension, namely as a Sine-Gordon model ( |1 . 1|) with the 
frequency j3 determined by the relation 

A £ ( P )=/3 2 /8tt. (5.4) 



Note that, since A f (0) = 1/2, the decoupling point corresponds to (3 — yAit. For this 
value of the coupling, the SG model is equivalent to a free Dirac fermion (see (|2.6|)), as 
expected from the well known equivalence between the thermal Ising model and a free 
Major ana fermion. 
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Consider now a generalisation of the Ashkin- Teller model defined by the action 

Ah = Aat + h / c/ 2 x <7i(x)<7 2 (x) , (5.5) 



with an additional spin-spin interaction between the two planar Ising models. 

For t = 0, it can be shown that the action (|5.5|) defines an integrable massive field 



theory whose spectrum and scattering amplitudes have been determined in |23 . 

For r^O, there is however the possibility to make a fine-tuning of the two parameters 
r and h in such a way that a massless RG flow is obtained. The infrared fixed point can 
be easily identified by means of the following physical consideration. Let us consider 
first the strong-coupling limit h — > oo: in this regime, the last term in ( |5.5| ) forces 
the spins <i\ and cr 2 to assume the same value on each site. Hence, the system results 
effectively reduced to a single Ising model. The reduced system is generally massive but, 
as any other Ising model, it may become critical provided the temperature of the two 
original layers is appropriately tuned to its critical value. It can be easily argued that 
the strong-coupling limit scenario just described should also occur for finite values of h: 
the coupling h tends to order the system (hence introducing a finite correlation length) 
but this tendency can be contrasted by a value of the temperature sufficiently high. In 
conclusion, for each value of p, there should be a critical line in the r/i-plane along which 
the system flows from the C = 1 fixed point to the C = 1/2 Ising fixed point. For the 
symmetries of the problem, the approach to the Ising fixed point must occur along one 
of the irrelevant directions of the conformal family of the identity operator: in fact, the 
massless flow is along the line which separates the high- and the low-temperature phases, 
hence the conformal families of the magnetic and energy operators of the Ising model are 
both ruled out since the former is odd under the Z 2 spin symmetry while the latter is 
odd under the high-low temperature duality of the model. So, only the conformal family 
of the identity operator is left. 

Let us finally rephrase the above discussion in the bosonic language. Since £ ~ cos /3<p, 
eq. ( |5.3| ) suggests the identification P ~ cos(f3ip/2 + 5p), with 5p to be fixed. The model 
( |5.5| ) then corresponds to the two-frequency SG model ( |3.1| ), with the identification 
\i = r, A = h, a — [3/2, 5 = 5p, with j3 determined in terms of p through ( |5.4|) . We found 
in the previous section that, for a/ (3 = 1/n, a flow to the Ising fixed point takes place in 
the A/i-plane if \5\ = ir/n. This implies \5p\ = tt/2, so that P ~ sin(/?y?/2). 
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5.2 Massive Schwinger model 

The massive Schwinger model describes quantum electrodynamics in two dimensions. It 
is defined by the Lagrangian 

Cmsm = -\f^F, u + ^Yd, -m- erfA^ . (5.6) 

If the "quark" mass m is set to zero, the theory reduces to the ordinary Schwinger model 
which is known to be equivalent to a free scalar field of mass e/y/7T p6| , On the other 
hand, if we set e = in ( |5.6| ) we are left with a theory of free massive fermions which, as 
already remarked in section 2, is equivalent to a pure SG model with f3 = \/A~k. Putting 
this information all together, one obtains the following bosonised form of the theory Q5.6D 



1 e 2 

C-Bose = o(<^) 2 ~ 7r(v + l) 2 + cm 2 cos VAttlp , (5.7) 

where c is a positive constant whose precise value is immaterial for our considerations. 
The parameter 7 accounts for the fact that, contrary to what happens in the two limiting 
cases m = and e = 0, a shift of the field p does not leave invariant the complete 
theory (|5.7|) : it is in fact proportional to the angle of the gauge lagrangian (|5.6j ) which 
physically corresponds to the presence of a constant background electric field which breaks 
parity |25| . 



The bosonic lagrangian ( |5.7|) can be seen as the specialisation of the theory (|3.1|) to 
the case (3 = a/4~7t, A = (e/v^ra) 2 , 5 = ccy, in the limit a — > 0. The analysis performed 
in the previous sections for the theory (|3.1|) is easily adapted to the present case. The 
requirement that the lagrangian ( |5.7|) describes theories which do not differ simply for 
a shift of the vacuum leads to the restriction I7I < a/vt/2. Chosing the G angle to take 
values in the interval (— 7r,7r), then one has 7 = Q/^/Att. The potential term of the 
lagrangian ( |5.7|) for a generic value of G and for the special value G = n is drawn in 
Figure 8. 

Let us apply the perturbative arguments of Section 4 to the limit e/m — > of ( |5.7| ). 
The unperturbed theory is the SG at the free fermion point. The spectrum consists of 
soliton and antisoliton only, which in the fermionic description must be identified with 
the non-interacting quark and antiquark. The perturbing operator \1/ = (9? + 7) 2 is non- 
local with respect to s and s. Such a non-locality is of a less simple form than that of the 
exponential operators e*"^, but is however easily deduced from the latter. In fact, taking 
the k-th derivative with respect to a of eqs. ( |A.4|) and ( |A.10| ) and then setting a = 0, one 
finds 

(O|/(O)| S (^ + 2z7r),s(0 2 )) = (0|(^(0)-27r//3) fe |^ 2 ), S (^)) , (5.8) 
-zRes 9l _ e2=±m (0|^(0)| S (^),^ 2 )) = (0|v? fc (0) " (<^(0) ^ 27r//3) fc |0> . (5.9) 
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With j3 = y/^n and (0|<^|0) = 0, the last equation gives 

iRes e= ± ln F*(9) =tit® ■ (5.10) 

The interpretation of this result follows from the general discussion of Section 4. For 
|0| < 7r the soliton and antisoliton masses receive an infinite correction. In the fermionic 
description this amounts to say that quarks and antiquarks are confined by an arbitrary 
small electromagnetic interaction. Quark-antiquark pairs are bounded by a linear poten- 
tial and give rise to a string of electrically neutral bound states. The confining potential 
and an estimate of the number of bound states as a function of energy are given by the 
formulae (Q and Qgg ) with A = e 2 /2 and B = Q/n. 

For |0| = 7r, however, the annihilation pole in the soliton-antisoliton form factor of 
the perturbing operator is (partially) cancelled and the mass correction is finite. Quark 
and antiquark survive unconfined in the perturbed theory as kinks interpolating between 
the vacuum located at (p = 0(e 2 ) and a degenerate adjacent minimum (the latter lies 
at ip = ±y/n + 0(e 2 ) depending on the sign of 6) (Figure 8.b). Since no topological 
excitation is present in the opposite limit e/m — > oo, there must exist a critical value of 
the ratio e/m for which a confining phase transition takes place. The arguments of the 
previous section apply basically unchanged to the present situation and one concludes 
that along the critical trajectory in the em-plane the system flows from the C = 1 to the 
C =1/2 fixed point. 

6 Conclusion 

The multi-frequency Sine-Gordon model is an interesting example of non-integrable 
QFT. In this paper we have exploited the Form Factor Perturbation Theory to study 
the evolution of the particle spectrum of this model as a function of the parameter rj 
which labels its interpolation between two integrable SG models. The soliton sector is 
the one which results the most affected by the non-integrable dynamics: as discussed 
in Section 4 and also illustrated by the examples in Section 5, its evolution strongly 
depends on the values assumed by the frequency ratio u> and the phase-shift 5. In the 
FFPT, disappearing of soliton states and confinement of their multi-particle excitations 
originates from the non-locality of the perturbation operator with respect to the field 
which creates the one-soliton state. In the special cases u — 1/n and 5 = ±l/n, a 
change of the vacuum state takes place by varying rj, with the consequent presence of 
phase transition lines: these are generally associated to massless field theories which 
interpolate between the conformal field theory with C = 1 and the one associated to the 
Ising model, with C = 1/2. 
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It would be highly interesting to match the rich scenario of the multi-frequency SG 
model as coming from the analysis of its particle spectrum with a more detailed inves- 
tigation of its dynamics, in particular the calculation of its scattering amplitudes. As 
shown in Ref. |28| for the double SG model, a remarkable pattern of resonances entering 
the collision events of the kinks already emerges at the classical level. A challenging open 



problem is therefore to see how the results of Ref. |28] can be translated at the quantum 
level and conversion amplitudes between different types of kinks, inelastic events and 
energy levels of the resonance states of the non-integrable model can be computed. 
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Appendix A 

At the special values of the coupling £ = 7r/n, n = 1, 2, . . ., soliton-antisoliton scattering 
in the Sine-Gordon model becomes reflectionless, and is described by the transmission 
amplitude 

Z-j sinh ~[0 + Z7r(l — k/ri)] 
S{6) = ~ 11 sinhlp-Mra-Jfc/n)] ' (A ' 1} 

in which the simple poles located at 9 = — k/n), k = 1, 2, . . . , n — 1, correspond to 
the breathers appearing as bound states in ss scattering. 

As discussed in Section 4, the soliton-antisoliton form factors of the exponential op- 
erators 

F%{6 X - 9 2 ) = (Ole^ls^)^)) (A.2) 
satisfy the unitarity and crossing equations 

= S(6)F«{-6) , (A.3) 

F«{9 + 2in) = e- 2lna/t3 F«(-6) . (A.4) 

Since we have chosen the vacuum of the theory to be located in ip = 0, the reflection 
ip — > —if maps s in s. The invariance of the theory under this reflection then requires 

KM = Ff«(d) . (A.5) 



We parametrise the matrix element ( A.2 ) as 



where 



satisfies 



Fq(9) = —i sinh - exp 



dxsinhf(l-e/vr)sin 2 fc^ 



2tt 



o x sinh f£ cosh § sinh a; 



(A.7) 



F o (0) = (~l) n+1 S(9)F (9) , (A.8) 

F (# + 2nr) = F o (0) . (A.9) 

The denominator sinhra^ in ( |A.6| ) contains the kinematical pole in 9 = m associated to 
the residue equation 

- iRes e=i ,F%{9) = (1 - e -W/3)( e ^> f (a.io) 
as well as the bound state poles associated to the residue equations 

- i Res e=in(1 „ k/n) F? § (9) = T k F« (A.ll) 
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(here r fc = [— i Res0 = j 7r (i_/ c / n )S'(6')] 1//2 is the three-particle coupling between soliton, anti- 
soliton and fc-th breather and F£ = (Ole 40 ^ - 1 !^)). The pole in 9 = is cancelled by the 
zero contained in Fq(6). 

According to eqs. (|A.3| ), ( |A.4|) and (|A.5|) , P a {9) entering (|A.6|) must satisfy 

P a {6) = (-l) n P- a (-9) , 
P a (9 + 2 i n) = (-l) n+1 P a (9) , 



and then can be written as 

'3= 
sM 



Pa(9) 



Ef=0 [%' 



a e 



jt 



-a e 



-iff 



n odd 
n even. 



(A.12) 
(A.13) 

(A.14) 



E-Lo M«) e °' +1/2)e + M-a)e~ ' +1/2) 1 > 
The integers N and M can be fixed from the asymptotic behaviour of the known solutions 
for a = ±/3 and a = ±(3/2 §. One finds N = {n - l)/2 and M = (n - 2)/2. It is 
convenient to define a (a) = ao( Q; ) ~ a o( — a )> Cj(a) = aj(a)/a (a), dj(a) = bj(a)/b (—a) 
and to rewrite P a {9) in the form 



P a (9) 



a {a) 



1 + Efc 1)/2 fe(a)^ + Ci (-a)e-^ 



6d(«) E5=o 2)/2 JSe (j+1/2)e + diC-aJe-^+W 



n odd 
n even 



(A.15) 



in which an inessential overall normalisation has been made explicit and (both for n odd 
and n even) n — 1 a-dependent coefficients appear as the only remaining unknowns in 
the problem. They are fixed exploiting the following observation. The lightest breather 
P>i, being the particle interpolated by the elementary field (p, changes sign under the 
reflection (p — > —if. The k-th breather can be seen as a bound state of k breathers B\, 
and then behaves as — > (— l) k B>k when if — > —if. As a consequence we have 

Fu = (-D k Fr a , 



k = 1,2 n- 1. 



Hence, remembering eq. QA.ll]) and defining Q a {9) = e a9 ^P a (9), we obtain 



Qc 



nk 



n 



(-l)*Q a -i 



1,2 n-1. 



(A.16) 



(A.17) 



This is a set of n — 1 linear equations which uniquely determine the coefficients in ( |A.15| ). 
The solution is 



7T 



cj{a) 



ntan — I m 
2n 

m=l 



n 



cot — 

n 



a 1 



71+1 



777. 



1. 



n — 1 



for n odd, and 

d i( a ) = n 

m=0 

for n even. 



7T 



2n 



77 



tan — m cot — 



7T 



n 



a 1 /n 
P + 2\2 



m 



0,1, 



n 



2 

(A.18) 
(A.19) 
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Appendix B 



We discussed in Section 4 how the topological excitations of SG model can undergo a 
confinement phenomenon under the action of a perturbation. In this Appendix we give 
a more quantitative description of the confinement pattern within the semiclassical ap- 
proximation. 

Consider at first a soliton-antisoliton (ss) state in the pure SG model fll.lj) . Classically, 
it is associated to a configuration of the field <p taking the vacuum value <p = from minus 
spatial infinity to a point xi(t) (the position of the soliton at time t) where it switches 
to the value <p = 2%/ '(5 of the right adjacent minimun of the potential; it keeps this 
value until the point x 2 (t) (location of the antisoliton) where it switches again to ip = 0. 
Similarly, the antisoliton-soliton (ss) configuration corresponds to an interpolation from 
if = to the left adjacent minimum ip = —27i/j3, and back (see Figure 9). 

Let us add to the SG potential the perturbing term V((p) = — Xcos(a(p + 5) which 
lifts the degeneracy among the vacuum at <p = and its two adjacent minima^]. The two 
configurations ss and ss described above are no longer stable since a surplus of energy 
equal, respectively, to U sS = (x 2 — xi)[V(2tc/ f3) — V(0)] and U Ss = (x\ — x 2 )[V(— 2n//3) — 
V(0)] is associated to the intermediate plateau. The tendency of the system is to shrink 
in order to minimise the energy of the configuration: hence, the new term gives rise to 
the following attractive potential between the soliton and the antisoliton 

U(x) = 6(x)U Ss + 9(-x)U s§ = A(\x\ - Bx) , (B.l) 

where 9(x) is the step function and 

7TCX ITOt 

A = 2\ sin — cos £ , S = cot — tan5 . (B.2) 

P P 

Such a linearly rising potential confines the soliton-antisoliton pair and gives rise to a dis- 
crete spectrum of bound states. A continuum spectrum corresponding to asymptotically 
free motion is recovered only in the limiting case \B\ = 1. 

The number of bound states N(E) as a function of energy can be obtained by the 
semi-classical formula 

N(E) = — (ppdx , (B.3) 
2ix J 

where the integral is computed along the classical orbit. In a non-relativistic approxima- 



17 For the reasons discussed in Section 3 the perturbation must be such that the central mininum 
remains the lowest among the three we are considering here. This amounts to require A positive and 



B\ < 1 in (B.l) below. Since the present considerations involve only a localised region of the potential, 



they lead to less severe restrictions on S than those obtained in Section 3. 
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tion, the above integral computed with p = [M(E — 2M — U)] 1 ^ 2 gives 

AVM(E-2M)^ 
N(E) ~ 3^(1 - W) ' (R4) 

where M is the soliton mass and the rest energy 2M of the soliton-antisoliton pair has 
been included. The number of stable bound states is N(Et), Et being the value of the 
lowest threshold in the energy plane (Figure 10). Et equals 4M for j3 2 > Att and twice 
the mass m\ of the lightest breather for f3 2 < Att. Comparing with the exact solution 
expressed in terms of Airy function of the Schrodinger equation with the potential ( [B.l| ), 
the semi-classical formula (|B.4|) provides excellent estimate of the energy levels (an error 
of 0.7% for the first energy level, which is further reduced at 0.1% for the second one). 
Therefore, for all practical purposes, we can consider formula (|B.4|) as an exact one. 



The computation of the integral ( |B.3| ) with p = [(E — U) 2 — M 2 /4] 1 / 2 required in the 
relativistic case cannot be performed exactly^. However, it can be checked numerically 
that the non- relativistic computation provides a reasonable estimate for N{Et)- 

The stable bound states originating from the threshold of soliton-antisoliton config- 
urations may influence the large asymptotic behaviour of correlation functions. This is 
certanly the case for those fields which were coupled in the unperturbed system to the ss 
state. Let 0(x) be one of these fields. Neglecting states with higher number of particles 
and masses, the two-point function of the operator 0(x) in the unperturbed system was 
given by 



2 



G{x) = (O(x)O(O)) ~ / ^-p(fi 2 )K (fi | x |) (B.5) 
with Kq(t) the Bessel function and 



(2M) 2 27T 



Hence, the leading behaviour for large | x \ in the unperturbed case is given by 

G(x)~D J-—^— (B.7) 

V 7T \ x 2 +2<T 

where D and a are defined by the behaviour of the spectral function p(/i 2 ) near the 
threshold 2M , p(/i 2 ) ~ D (/i — 2M) 2a . When the branch cut associated to the threshold 
of the state ss breaks down into the infinite sequence of bound states - iV& of which stable 
and whose energies are obtained by inverting eq. ( |B.4| ) - the correlation function becomes 

4 Nb 1 

n=l 
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It is easy to see though that in this case, N(E) ~ E 2 for large energy. 
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Hence, the leading asymptotic behaviour for | x |— > oo is determined by the lowest of 
them and we have 
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Figure Caption 



Figure 1 . Topological excitations: (a) soliton state, (b) anti-soliton state, (c) multi- 
soliton state. 

Figure 2 . (a) Perturbation of the SG potential (full curve) by another periodic in- 
teraction (dotted curve); (b) when | 8 |> the vacuum moves away from the 
origin. 

Figure 3 . Examples of degeneracy of the vacuum at the origin with one of the other 
{n — 1) minima for 5 = -: (a) uj = — (m ^ 1); (b) u = -. 



Figure 4 . Straight line ( 3.1 5| ) on the torus T in the rational (a) and irrational case 



(b) . The black point identifies the location of the minimum of the potential on the 
torus. 

Figure 5 . Form of the potential near the phase transition point: (a) i] < i] c ] (b) i] = r] c ; 

(c) 7] > 7] c . 

Figure 6 . RG flows from C = 1 to C = \. 
Figure 7 . The two-layer Ising model. 

Figure 8 . Potential term in the massive Schwinger model: (a) ^ ir; (b) = 7T. 

Figure 9 . (a) Two-particle kink-antikink and (b) antikink-kink states. 

Figure 10 . Bound states of the linear potential confining the soliton-antisoliton pairs. 
The levels above E? correspond to unstable particles. 
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